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Abstract 

The non-relativistic bosonic ground state is studied for quantum N- 
body systems with Coulomb interactions, modehng atoms or ions made 
of N "bosonic point electrons" bound to an atomic point nucleus of Z 
absolute "electron" charges, treated in Born-Oppenheimer approxima- 
tion (the nuclear mass M = oo). By adapting an argument of Hogreve, 
it is shown that the (negative) bosonic ground state energy 8,''^{Z,N) 
yields the monotone non-decreasing function N i— )• E.'^{XN , N) / for 
any A > 0. The main part of the paper furnishes a proof that when- 
ever A > A, ?s 1/1.21, then the limit e(A) := limN ^oo S.-^{XN, N)/N^ 
is governed by Hartree theory, and the rescaled bosonic ground state 
wave function factors into an infinite product of identical one-body wave 
functions determined by the Hartree equation. The proof resembles the 
construction of the thermodynamic mean-field limit of the classical en- 
sembles with thermodynamically unstable interactions, except that here 
the ensemble is Born's, with as ensemble probability density func- 
tion on R^^, with the Fisher information functional in the variational 
principle for Born's ensemble playing the role of the negative of the 
Gibbs entropy functional in the free-energy variational principle for the 
classical petit-canonical configurational ensemble. 
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1 Introduction 

In Spring 2011 Brookhaven National Laboratory announced the discovery of 
the first dozen and a half anti-a particle^] [STARllj . Since a particles in 
their ground state have spin zero, the same can be expected for anti-a parti- 
cles which, therefore, will be twice negatively charged bosons for practically 
all quantum-mechanical low-energy phenomena in which they participate. In 
particular, it doesn't take much imagination now to predict that lab-produced 
atoms with such bosonic "electrons" in place of the usual fermionic electrons 
are eventually going to become an experimental reality, although the path 
to their controlled production and storage is cluttered with enormous tech- 
nical obstacles that need to be overcome. Yet such a feat should be possi- 
ble, as demonstrated by the recent success story about storing anti-if atoms 
|HHetall2] . and by the recent lab creation of "protonium" |Zetal06al IZetalOGb] . 
an atom made of a proton and an anti-proton, both fermions. 

The simplest atom formed with an anti-a particle as electron substitute, 
and presumably the first one to be produced in the lab, would seem to be 
"alphium," i.e. the a particle analog of protonium and of the familiar positro- 
nium (which I'd much rather like to see called "electronium" ) . Alas, although 
an (anti-a,a) bound state is the simplest special case of an atom with bosonic 
"electrons," a bosonic atom having just = 1 bosonic "electrons" cannot dis- 
play any of the effects which depend on the permutation symmetry of bosonic 
many-body wave functions. To see those, laboratory-produced atoms with 
a-bosonic "electrons" would have to be formed with A^ > 1 bosonic anti-a 
particles attracted by a conventional nucleus of charge 2eA^. So the atom with 
N = 2 a-bosonic "electrons" would be the simplest truly bosonic atom; it'll 
have a Beryllium nucleus, the stable isotope of which (^Be) is a fermion with 
spin 3/2. In principle it is conceivable to go up to A^ ~ 46, requiring a Ura- 
nium nucleus, although it would seem exceedingly difficult to strip away all 
electrons of the atomic hull surrounding any heavier nucleus. 

Of course, alphium and the heavier bosonic atoms with anti-a particles 
as electron substitute will be short-lived, their anti-a particles annihilating 
with (parts of) their nuclei in a complicated fashion, thus opening an inter- 
esting new venue for studying the strong interactions at their lower energies. 
However, if positronium and protonium are any indicators, then alphium and 
the heavier a-bosonic atoms may well have lifespans of the order of //-seconds. 
This is long enough to gather empirical data about electromagnetic transitions 
between their excited atomic states and their atomic ground state, for which 
the treatment of anti-a particles as bosonic "electrons" is vindicated. 

^This announcement came very much to the delight of the author, who in 2009 jKie09b] 
wrote: "Since fermionic anti-'^He nuclei have already been produced in heavy ion collisions 
at CERN [Aetal96[ IAetal03al lAetalOSb] , it seems a safe bet to predict that also bosonic 
anti-a particles are going to be produced in the laboratory,...", but who had no idea that 
this was such a timely assessment of the experimental situation. 
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Data on electromagnetic transitions in these a-bosonic atoms should not 
only yield information which depends on the symmetry of many-body bosonic 
wave functions, but also information about the structure of the anti-a particles 
themselves, for these are not at all point-like, in contrast to the conventional 
fermionic electrons. For instance, the fine structure of the a-bosonic atomic 
spectrum should show some novel effects due to electric quadrupole-quadrupole 
coupling, which is absent in atoms with true electrons. 

No later than when atoms with anti-a particles as bosonic "electrons" are 
being produced in the laboratory will theoretical physicists get busy calculating 
the various annihilation channels (a challenging task) and the electromagnetic 
spectra of these bosonic atoms (a more "straightforward" task). The calcu- 
lation of the bosonic spectra will require accurate numerical approximation 
methods, but very likely also the traditional Hartree approximation combined 
with perturbation theory for Schrodinger operators. Happy news for mathe- 
matical physicists who have been busy with rigorous studies of such atoms for 
some time already; see |BeLi83j . pe&O] . [SoIOO] . |Bac91j . |BLLS93j . |Rus97j . 
[BaSeOlj . |Kie09b] . |Hogll| , and in particular |LiSelO] and references therein. 

In this work we will concern ourselves with the Hartree approximation to 
the ground state of a non-relativistic atom or ion with bosonic "electrons." 

Informally, we can summarize our results as follows: let E,'^{Z,N) with 
Z = XN denote the bosonic ground state energy of a Born-Oppenheimer atom 
(A = 1) or ion (A 7^ 1); the index "00" stands for the infinitely massive nucleus. 
Let tjji^^ denote the (suitably scaled) bosonic ground state wave function of 
this iV-body system, i^(0) the pertinent asymptotic (iV- independent) Hartree 
functional of a one-point wave function, and A^'' the normalized empirical n- 
particle measure (as obtained in an ideal measurement of the positions). 
We will show that whenever A > A* ~ 1/1.21, then: 
(i) N-^E%{\N, N) / inf^ :^a(0) when ^ 00; 
{ii) ^x{(f)) has a unique minimizer, (px, modulo sign; 

(2m) ^i^^ ^ 0^ as AT ^ 00; 

(iv) Aj^^ |0Ap ® ■ ■ ■ ® |0Ap (n-fold product), as N 00. 
The precise statements are given in section II, their proofs in section HI. The 
remainder of this introduction relates our work to the literature on the subject. 

First of all, {{), (ii), {in), and (iv) only express what everyone should expect 
to be true anyhow; in fact, (z) and [ii) are not even new. Namely, except for the 
monotonic increase, (i) has been proved in |BeLi83] : the monotonic increase 
was proved, for A = 1, in |Hogll| . As for (ii), this was shown for A > A* in 



|BBL81j : the numerical value A* f« 1/1.21 is due to |Bau84j and |Sol90j . Yet, 
our proof of (i) and (ii) overall is novel (to the best of the authors knowledge). 

As to (iii), this result seems to be new, though superficially it is reminiscent 
of the results proved recently in |BGM00j . [ErYaOlj (see also jBEGMY02j ). 
[ElSch06j ■ |FKS09] . and in |Picll] . which establish the asymptotic exactness 



3 



of the dynamical Hartree equations for Coulomb or Newton interactions in the 
limit N ^ oo; see also |Hep74] and |Spo80] for earlier results with more regular 
interactions. In these works it is shown that ipl^^ — > (j)f^ as — ?■ oo z/ tpQ^^^ — )■ 
0®^ as — )■ oo, where ipl-^'' denotes the A^-body Schrodinger wave function at 
time t and ^/'g^'* its initial state, while (pt is the solution to the Hartree equation 
at time t launched by initial data (pQ. These "propagation of factorization" 
results are very much in the spirit of what Kac, McKean and Lanford called 
"propagation of chaos" in their approaches to derive Boltzmann's equation in 
the kinetic theory of a dilute classical gas. Note that the factorization of the 
initial wave function (as — )■ oo) is being assumed in these works; in fact, 
most A^-body wave functions do not factorize in the limit N ^ oo. 

Also (iv) is new. Previously, in |BeLi83] and |LiYa87] (taking the c — t- oo 
corollary) , it was shown that J \ilj{^^\'^{q^, ...,qj^)d^q2 - ■ -d^qi^ |0a(?i)P- Our 

{in) generalizes this to J \ilj[''''\'^{q^, ...,q^)d^q^^^ - ■ ■ d^q^ Ilk=i\<P^i(lk)? 
for any n G N. Such results state that the expected value of the normalized 
empirical ra-particle density agrees with (|(/>aP)®"' in the limit A^ — )■ oo. The gap 
between such an average result and (iv) is bridged by a law of large numbers. 

Our last statement makes it plain that our approach will be probabilistic 
(or: statistical-mechanical) in nature. More precisely, it resembles the con- 
struction of the thermodynamic mean-field limit of the classical ensembles 
with thermodynamically unstable interactions, except that here the ensem- 
ble is Born's, with IV^a'^^P as ensemble probability density function on M^^. 
The proof exploits the fact that the finite- A^ variational principle for I^^^^P 
resembles the free-energy variational principle for the classical petit-canonical 
configurational measure, with the Fisher functional in the variational prin- 
ciple for Born's ensemble now playing the role of the negative of the Gibbs 
entropy in the classical petit-canonical configurational free-energy. Readers 
familiar with the works MeSp82] , |Kie09aj . |Kiellj on the petit- and micro- 
canonical classical ensembles will recognize the strategy; see also |APT94j for 
a petit-canonical quantum ensemble which uses the classical Gibbs entropy 
for distributions on Feynman-Kac path configuration space. 

We end this introduction on a cautionary note. The ratio of the a-bosonic 
"electron's" mass m versus the mass M of any regular nucleus is not really so 
small that one could have much confidence in the Born-Oppenheimer approx- 
imation (M — 7- oo). The correct Galilei-invariant atomic (or ionic) model with 
a dynamical nucleus of mass M is a more subtle A^ + 1 body problem. Slightly 
simpler than this is the Galilei-invariant Newtonian gravitational analog of the 
problem, [Pos62j , [Hal92j . Its pseudo-relativistic version, known as a "bosonic 
neutron star," has recently attracted lots of attention, see [LiYa87j . |ElSch06j . 
[HaLuOGj . IFrLeOTj . |HaLu08j . |FrLe09] . |FKS09] . and [PicTT] . We hope to 
address such Galilei-invariant systems in a follow-up work. 
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2 Statement of results 

In this section we give a precise formulation of the A^-body variational prob- 
lem (subsection 2.1) and of our Hartree limit theorems (subsections 2.2 and 
2.3). The bosonic ions are treated in Born-Oppenheimer approximation, for 
simplicity. Also, all particles are treated as pointlike, which can only be a 
rough approximation for a-bosonic atoms and ions. 



2.1 The A/^-body variational principle 

The Hamiltonian of a non-relativistic atomic ion with an infinitely massive 
nucleus (the Born-Oppenheimer approximation; indicated by the label oo) is 
the formal Schrodinger operator 

Ht'''= E (^IP.P-^-^e^^)+EE^V^, (1) 

l<fc<iV ^ ^ l<j<k<N ^ 

with the understanding that J2J2i<j<k<i ~ ^- ©> Pk = —iKVk is the 
familiar momentum operator canonically dual to the A;-th "electron's" con- 
figuration space position operator, q^. e R^. Moreover, m is the "Newtonian 
inertial mass" and ze the electric charge of each of the N & N "electrons," with 
z = —2 for anti-a particles and z = —1 for true fermionic electrons. Lastly, 
Z\z\e (with \z\Z e N) is the charge of the atomic nucleus (which is fixed at 
the origin); of course, e(> 0) is the conventional elementary electric charge. 
Although I^IZ G N in nature, we will conveniently allow Z G M+. 

The formal operator h!^'^^ is densely defined on (£[f fl £^)(M^^). As 
self-adjoint extension we take its Friedrichs extension |ReSi80] . also denoted 
by H^'^\ which is a permutation-symmetric, self-adjoint operator with form 
domain given by the A^-fold tensor product D^^^^ = S)^(R'^) ® ■ • • ® iD^(]R^) C 
11^ (M'^^). The quadratic form associated to the Friedrichs extension h!^'^^ is 

^,M)^^(N)^ ^ _ Zz^e^ c(^)(^('^)) + z^e^ 7(^)(^(^)), (2) 

where 

^(iv)(^(iv))^ r ^ |v,vi/(-)|2d'^, (3) 

J l<k<N 



J l<k<N 



(4) 

J l<k<N 

/(-)(v^(-))= [ j: j: j^^i^'^'W- (5) 

J l<k<l<N 

The bosonic ground state energy of H^x,'^^ is defined by 

£^(Z,Ar):=inf{Q(^'-)(v[/(-))|v[;(-)GD^"^ ||v[/(-)||^.(^3.) = l} . (6) 
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The infinimum is always finite, yet depending on the ratio of Z vs. it may 
or may not be achieved by a minimizing A^-body wave function ^%'''^ G D^l\ 
called the ground state. Thus, the "ground state energy" may or may not be 
the "energy of a ground state." 

To sort this out we recall that ^'^{Z, N) < min aess (-ff^'^'*), by the min-max 
principle (Theorem XIII. 1 in |ReSi78] ). So either £^(^, N) = min cxess 
or E-g{Z,N)<mmaess{Ht''^). 

When E,^{Z,N) = min cTess (-f^oo '^'*) , then a minimizer ^I/l^'^"* may or may 
not exist, indeed (to the best of the authors knowledge). However, the bottom 
of the essential spectrum of a general self-adjoint operator can only be one 
or more of the following: (a) an eigenvalue of infinite multiplicity (e.g. the 
eigenvalue of a projection operator onto a finite-dimensional subspace of 
Hilbert space), (b) a limit point of the pure point spectrum (as is the case for 
-f^oo'^''), or (c) the bottom of the continuous spectrum (as is the case for Ht'^ 
also; note that the possibilities are not mutually exclusive). Clearly, in case (a) 
the ground state would be infinitely degenerate, and if a ground state exists in 
cases (b) or (c), even if different from ), it would be "as close as it gets" 

to being infinitely degenerate without literally being infinitely degenerate; we 
will call this infinitely quasi-degenerate. Paraphrasing Simon |Sim70] . all such 
ground states "are somewhat pathological beasts." 

On the other hand, when E''^{Z,N) < min(Jess(-f^oo'^''), the existence of 
a minimizing ^E'^'^'' is guaranteed by the min-max principle. In this case 
is in fact nondegenerate and it can be taken to be real-valued and 
positive. Moreover, the permutation symmetry of h!^'^^ implies that '^^^''^^ is 
automatically permutation symmetric, too, i.e. bosonic. All this is of course 
well-known; see, e.g., |Lie90] . |LiSelO] . 

For our approach we will need the existence of a nondegenerate '^^^''^^ G 
D^i^\ So we follow |Sim70] and consider to be a proper ground state if 

and only if ^'^{Z, N) belongs to the discrete spectrum, equivalently E,''^{Z, N) < 
min (Tess {H^''^^) . Henceforth ^l/^'^"* always denotes such a proper ground state. 

The condition E'^{Z,N) < min (Tess (-f^oo''^'') can be expressed in a more 
user-friendly manner. This is accomplished by the HVZ theorem, which asserts 
that aess(^^''^^) = [E,"^{Z, N -l),oo), see Theorem 11.2 in [TesOQ] (this is a 
special case of the HVZ theorem for H^''^'' with nuclear mass M < oo, see 
Theorem XIII. 17 in |ReSi78] ) . So a proper ground state ^l/^'^"* exists whenever 
£;^(Z, A^) < 8.'^{Z, N — 1). We need to know for which combinations of Z and 
this inequality is fulfilled. 
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2.1.1 The nuclear charge Z necessary to bind bosonic electrons 

By Zhislin's theorem (Theorem 12.2 in jLiSelOj ). a proper A^-body ground 
state ^l/^'^'* exists whenever Z > N — 1. On the other hand, Benguria and 
Lieb |BeLi83] have shown that exists whenever Z > X^N provided N is 

large enough, with = 1/(1 + 7) iii the notation of |BeLi83] . and numerically 
7 ^ 0.21 |Bau84] . Combining these two results we obtain the following. 

Corollary 1. There exists a X* satisfying 1 > A* > A* such that for all N & N 
a proper ground state wave function \E'^^''^' for ^ exists whenever X > X* , 
and there exists some N* G N such that for all N > N* a proper ground state 
wave function \t'l^'^''^' for ^ exists whenever A > A* ~ 1/1.21. 

Remark 1. In view of CorollaryUlwe will have to restrict N to be "big enough" 
in the proofs of our main theorems, which are stated for the optimal parameter 
region A > A*. Whenever X > X* , this restriction on N can be dropped; in 
particular, this is the case if X> 1 . 

We emphasize that any nontrivial values of A* and A^*, i.e. values other 
than A* = A* and A^* = 1, whatever they might be, will presumably just be 
artifacts of our incomplete knowledge about the full range of Z values, given 
A^, for which a proper bosonic ground state exists. Indeed, as made plainly 
clear by the results of |BeLi83] . Zhislin's theorem does not reveal the full range 
of such Z values. We next state a conjecture as to what this range could be. 

We define Z^{N), the minimum amount of nuclear charge needed (in the 
sense of: to be exceeded) to properly bind A^ "bosonic electrons," thus: 

Z^iN) := mi{Ze M+IS^g'^^G^lf ^s.t (^'^H^S'"^) < ^^iZ,N-l)}. (7) 

Remark 2. A ground state may or may not exist when Z = Z^,{N); however, 
such a ground state would belong to the essential spectrum and, therefore, not 
be proper in the sense stipulated above. Incidentally, it is even conceivable that 
non-proper ground states exist for Z < Z^{N). 

By Zhislin's theorem, and the virial theorem, we have < Z^,{N) < N — 1 
for all A^ G N; in particular, ^^,(1) = (the hydrogenic problem), and Z^{2) < 
1 (the Helium- type problem; discussed extensively in section 4.3 of |Thi02j ). 
Also recall that Zhislin's theorem states that a proper ground state '^^^''^^ 
exists whenever Z > N — 1. Thus it is suggestive to suspect the following. 

Conjecture 1. A proper N-body ground state wave function "^^^'^^ for ^ 
exists whenever Z > Z^{N). 



Remark 3. ConjectureUlis the analogue of a claim in ^ReSiTSJ, made (without 



proof) after their proof of Theorem XII. 9 in their discussion of "example 1 

r(o,N) 
'00 



revisited. " Indeed, just replace their Hq (= —A) by our He, 
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We note that the function i-> Z^:{N) is a cousin of Lieb's and Benguria's 
Z A^jjja^x(Z), which is the maximum number of (here: bosonic) electrons 
which a nucleus of charge Z can bind. Interestingly enough, the precise ana- 
logue of Conjecture [1] for Z i— )■ A^max(2^) is also an open problem. 

Conjecture 2. A proper N-body ground state wave function for ^ 

exists for all 1 < N < N^g^^{Z). 

Remark 4. Conjecture [H is phrased in ILiSelO^ more generally thus: If a 
nucleus with Z charges binds N electrons (and N > 2), then the same nucleus 
also binds A^— 1 electrons. This seemingly obvious "truth" has only been proven 
to be true for N < Z + 2, by Zhislin's theorem, for both bosonic and fermionic 
electrons. I thank the referee for drawing my attention to this conjecture, which 
is mentioned as an open problem in Chpt.l2 of \LiSelO^ . 

To state our last conjecture we define A,,(A^) := Z^.{N)/N. It follows 
from the work in |BeLi83j and |Sol90] that limjvtoo A*(A^) = A, ^ 1/1.21. 
Moreover, A*(l) = and < A^,(2) < 1/2. The upper estimate on A*(2) 
follows from Zhislin's theorem, which yields the monotone increasing upper 
bound A,,(A^) < 1 — for all e N. We suspect that A^,(A^) itself is monotone. 

Conjecture 3. The map N ^ ^*{N) is monotonic increasing. 

Remark 5. Conjectures\^ and\^ if true, would imply that N* = 1, so that in 
the proofs of our main theorems we could drop the restriction N > N* . 

We now move on to state our Hartree limit theorems. 
2.2 The Hartree approximation 

The Hartree approximation to the bosonic ground state energy 8.'^{Z,N) is 
obtained by estimating Q^c^'^\'^^^''^'') from above with the help of convenient 
trial wave functions = with $ G ^^(M^) satisfying ||$||£2(k3) = 1. 

We thus have £f (Z, A^) < inf$ Q,^''^^($®^) = inf$ :K^''^^($), where 
ry^(z,M) _ 5^(1) _ p^zz^e^ C^i) ($) + N{N - l)z^e'\ Z^^) ($ ® $) (8) 

is the Hartree functional for H^''^\ 

We are interested in the large- A^ behavior, with scaling Z = XN. Setting 
$(q) = N^/^(t){Nq) in (E]) yields J{^'^''^^(<l>) = N^^^{^^x{(t)) + O^N"^), where 

^,a(0) = ii'K^'\<P) - C«(0) + z\^\ 7(2) (0 ® 0) (9) 

is the "asymptotic Hartree functional," and 0{N'^) = — A^^z^e^^ /(^^(^ ® 0) 
a negative correction term. Dropping this negative correction term results in 
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the uniform (in A^) upper bound ^E-^{XN, N) < inf,^ :?4d,a(0) V A^, where 
the infimum is taken over G satisfying ||$||£2(]g3) = 1. Although 

this bound on N-^8,-^{XN, N) is weaker than the bound A^"^£f (AA^, A^) < 
Ar-3inf$J{(^^'^)($), the omission of the correction term worsens this upper 
bound only by an amount of relative order A^~^, which vanishes as N ^ oo. 

The Hartree approximation to the bosonic ground state energy, indeed to 
the bosonic ground state itself, becomes exact as A^ — > oo, provided A > A* : = 
limA*(A^). This is made precise by the following theorem. 

Theorem 1. (Asymptotic exactness of Hartree theory for bosonic ions.) 
Whenever A > A*, then 

{{} Ar-3£f (AAT, N) / inf : G & ||0||£2 = 1} as at ^ oo. 

Moreover, 

(a) 54o,a(0) has a unique (positive) minimizer, (px ^ -^^ H {||0||£2 = 1}, 
and then, setting il)'^^\q^, ...,q^) := N~i^'^^^^'^\N~^q'^, ...,N~^q^), we have 

(zz.) V^f ) <PT, 
in the sense that, for any n G N, weakly in 2} fl £3n/(3n 2) have, 

(m)'/|^?)pte,...,g^)d3g„^i---d3g^"-:iT R \Uqu)?- 



Remark 6. Except for the monotonicity, item (i) has previously been proved 



m 



' BeLiS^ . and so has {Hi)' in the special case n = 1; and item {ii) has previ- 
ously been proved in fBBLSlf . As stated for all n eN, item (Hi), viz. {iii)' , is 
novel. Also, the proof of TheoremUl overall seems to be novel. It is an "infor- 
mation variation on the entropy theme played" in lMeSp82^ , lCLMP92/9^ , 
IEM^, IAPT941 IChKiOO^ . lKie09a!^ . IMeTl^ . IKiWalS^ . 

Remark 7. The monotonic increase of N N~^E.^{\N, N) is a novel re- 
sult, too, strictly speaking; however, Hogreve \HoglVl already proved for neutral 
Born-Oppenheimer atoms that 8.^{N,N)/N^ grows monotonically in N >1, 
and inspection of Hogreve 's proof, see his formulas (7)-(10), reveals that it 
generalizes verbatim to Born-Oppenheimer ions, without restrictions on A. 

A slightly weaker atomic result, that 8,^{N, N)/N'^{N — 1) grows mono- 
tonically in N > 2, was earlier proved in IKieOQbf . Also this result general- 
izes to ions, with proper restrictions on A in place, i.e. for A > A* the map 
N I—)- 8,^{XN, N)/N'^{N — 1) is monotonic increasing. Of course we are hap- 
pier with Hogreve's result; also, his proof, which uses the bosonic symmetry of 
the wave functions, is simpler than the one in \Kie09bf . On the other hand, the 
proof in IKieOQb^ . which does not use the bosonic symmetry of the ground state 
wave functions, implies also that E.'^{N , N) / N"^ {N — 1) grows monotonically 
in N >2, where 



8,Z{N,N) :=min{Q(^)(vi/(-)) 



£2 



], (10) 
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with Dq^ = S^^lM."^) A • ■ ■ A i3"^(]R^), is the fermionic ground state energy of 

H^'^\ Since the fermionic ground state for / fiOj) is anti- symmetric by 
construction, it is not clear whether 8,'^{N, N)/N^ grows monotonically in N . 

2.3 Born's statistical ensemble 

To state our next theorem, we recall Born's statistical interpretation of Schro- 
dinger wave functions. According to Born, \il)'^^^\'^{q^, ...,q^) is the joint A^- 
body probability density for finding particle 1 at q^, particle 2 at ^^^d so 
on, when one performs an ideal measurement that yields the positions of the 
"electrons" in state ^Z^^^^. More precisely, let Q^^ be the random position which 
shows up in an ideal position measurement of particle /c, let Q*-^-* G be the 
random vector whose fc-th component triple is Q^^ and similarly let q^^^ G M?^ 
be the vector value taken by Q*-^-*, then for the ion's bosonic ground state, 

Prob (g(-)Gd'g near g(-)) = |V^f)|2(g„...,g^) n (11) 

^ ' l<k<N 

For each n<N, any Q^'^^ maps uniquely into the empirical random ?7-statistics 

A^QU^i,--;sn) = cr' E---E n sqS'^), (12) 

l<fcl<---<fe„<]V l<j<n J 

where the sj G are generic points in physical space M'^; see |EiSch02] . If one 
factors out the permutation group Sn, then the mapping between Q^^^ and any 
of its empirical [/-statistics is actually one-to-one. These empirical [/-statistics 
are themselves random measures (normalized counting measures of random 
point configurations), and the first few of them have practical significance. In 
particular, A^|lr) (s) is the normalized "electron" density of the bosonic atom 
found in an ideal density measurement, or rather as computed from an ideal 
measurement of all positions of the bosonic "electrons." For finite N the 
empirical electron density will fluctuate from any one such ideal measurement 
to the next, but fluctuations subside when N ^ oo. More precisely, the 
following law-of-large-numbers-tjpe result is proved in section 3.2: 

Theorem 2. (Asymptotic absence of density fluctuations.) 
For A > A*, and for any n eN, when N — t- oo, then 

(iv) Aq^^n) |0Ap ® ■ ■ ■ ® |0Ap {n-fold product), 
in the sense that when "dist" means Kantorovich-Rubinstein distance, then 

(iv)' Prob (^dist (^^l), (I^aH®'") > e) ^ OVe > 0. 

Remark 8. The case n = 1 in Theorem\^ supplies the proper meaning for the 
notion c»/A^|0aP as (continuum approximation of) the "particle density." 
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3 Proofs of the Theorems 



3.1 Proof of Theorem [T] 

Our strategy is to adapt the construction of the thermodynamic mean-field 
limit of the classical petit-canonical configurational ensemble with thermo- 
dynamically unstable interaction^ to "Born's ground state ensemble." The 
proofs exploit the fact that the finite- variational principle for Born's en- 
semble resembles the free-energy variational principle for the classical petit- 
canonical configurational measure, with the Fisher functional in the variational 
principle for Born's ensemble measure now playing the role of the negative of 
the Gibbs entropy in the free-energy functional for the classical petit-canonical 
configurational ensemble measure. 

We give the proofs for bosonic Born-Oppenheimer ions in general. The 
neutral atoms with Z = N are included as special case. 

3.1.1 Ensemble reformulation of Theorem [1] 

The proof of Theorem [T] uses Born's statistical ensemble interpretation of the 
rescaled ground state wave function tp^^^ of the Hamiltonian ([T]) only to the 
extent that J d^q = 1, so that formally the permutation-symmetric 

can be "called" an ensemble probability density. Truly probabilistic 
arguments will be employed eventually in the proof of Theorem [21 but the 
control for Theorem[T]is accomplished using functional-analytical and measure- 
theoretical arguments, only worded in statistical mechanics lingo. 

We write IV^J^^^P =: g['''^ e H S.^){R^^), by which we denote the 
permutation-symmetric probability measures on which are absolutely 
continuous w.r.t. Lebesgue measure. More generally, any normalized and 
rescaled bosonic iV-body wave function ip^^^ defines a formal probability mea- 
sure G n through j^/'^^^P =: g^''^ Now recalling that 
/ |V\E'pdg = JlVl^'lpdg, we see that upon defining Vs^'^Kli^ ■■■^ In) = 

we can express the quadratic energy func- 
tional (|2]) as a functional of q^'^\ Also pulling out a factor A^^, we get 
Q,(^-'-)(vl/{-)) = N''Ej^l{g(-)), with 

= i^''"') + -'eVf ) , (13) 

^See |MeSp82| for the original paper, where the classical petit-canonical configurational 
ensemble is treated for particles with Lipschitz-continuous interactions; see also [KieOQaj . 
[Kiell] for the classical micro-canonical ensemble with singular interactions. 

•^For optical reasons we prefer the notation over \/ {■)■ 
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where, for ^('^^ e ^^"(M^^^) with ^^^^^ e f|i(M3A'), 

=4 /|Vv/f?("f d'g (14) 



and 

vr = /(- AE i H- ^ E E E^)^'">d'?. (15) 

J ^ l<fc<]V l<k<l<N ' 

The functional 'Ej^^l^'*'^'') "^^^^ "^^^ quantum energy of It inherits 

from C^^'^^ the following properties: if A > A* and > A^*, the functional 
f fT3|) achieves its infimum at the unique and permutation-symmetric probability 
density q~^^ = i.e. 

if A>A* G [A„ 1), then holds ViV>l. Moreover, 

<i(^^?^)=iV-^ef(AiV,Ar). (17) 

We next draw the parallel to classical statistical mechanics. Namely, while 
f fTSj) equals the potential energy of g^^^ (up to the factor z^e^), with 'P' as usual 
standing for "Volta,"@ and f|T^ equals the kinetic energy of g^^^ (up to the 
factor hi^/Sm, with the strange factor 1/8 compensated for by the factor 4 in 
dH])), the rationale for the notation ^ is that f lT^ is identical to 

^(iv)(^(iv)) = j |Vln^(-)|%(-)d'g, (18) 

which was introduced by R. Fisher |Fis25] as measure for the information 
content of g^^^ regarding its mean. This Fisher information of g^^^ shares 
many properties with (the negative of) the Gibbs entropy of g^ \ defined bjl^ 

= /" In ((Vmc)=^^f?(^)) ^(^M'g (19) 



''We might have alternatively chosen C, to stand for "Coulomb," but we already used C 
to denote the expected value of the "central Schrodinger potential." 

^Normally in classical statistical mechanics the Gibbs entropy of a probability density 
function on phase space M.^^ x is taken relative to the uniform density h"^^ , which 
is the only place where one can make (a minimal) contact with quantum physics. Since 
we here work on configuration space M."^^ we do the next best thing and define the Gibbs 
entropy relative to the uniform density (h/mc)^^^ , where H/mc is the Compton wavelength 
of the "bosonic electrons." 
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if ^(^Mn^(^) e £^(M3^); otherwise is undefined (it could be ±00). 

Note that replacing by —5 in ( !T3|) gives us a formal classical configurational 
free-energy functional, with temperature mc^/'& (when lengths are normal- 
ized to h/mc), with c the speed of light in vacuo. Although that free-energy 
functional has no infimum, for the attractive Coulomb singularity of the cen- 
tral Schrodinger potential is catastrophic in three dimensionsj^ the analogy is 
compelling enough to suggest an adaptation of the strategy for constructing 
the thermodynamic mean-field limit of the classical petit-canonical configu- 
rational ensemble with thermodynamically unstable but bounded interactions 
|MeSp82] to "Born's ground state ensemble" for true Coulomb interactions. 
Indeed, we will now characterize the limit points of the sequence {Q\'^}n>n* 
in the spirit of the approach pioneered in |MeSp82] . 
Theorem 3. (Hartree limit of Born's ground state ensemble) 
For any ri G N and A > A*, let g^-^^ denote the n-th marginal probability 
measure of g^^\ Then the following holds: 

(a) weakly in (^3" fl £p)(M3") with p = 3n/{3n-2), we have 

lim Qi",! = pT; (20) 

(b) weakly in Sj^lM.^"-), we have 

lim VqZI = VpT-^ (21) 

(c) px G ^(M^) n {^/p G i3^(M^)} is the unique minimizer of the functional 

CooAp) ■■= £ J^'Hp) - ^VAC«(v/p) + z'e'l I^\^p ® >)); (22) 

(d) moreover, 

eoo,A(pA) = £oo(A) VA > A*, (23) 

where 

£oo(A) := lim iV-3£f(AiV,iV) (24) 

is a negative, strictly monotone decreasing, concave function of \> 

Remark 9. The functional Coo,a(p) defined in [2^) is just a rewriting of the 
asymptotic Hartree functional !Hx>,x{4>)j defined in in terms of y/p = (p. 
Therefore its unique minimizer p^ = \(f)x\'^, where (p^ minimizes the asymptotic 
Hartree functional. 



^Interestingly, the classical Coulombic free-energy functional is non-catastrophic in two 
dimensions, provided the temperature is higher than a critical value; cf. [CLMP92/95| , 
[Kii93] . [KlilT] ) 
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3.1.2 Proof of Theorem [3] 

We begin by stating an auxihary result which proves item (i) in Theorem [H 
Proposition 1. The map N h-)- infg(]v) j^'£!^\{Q'"^^) is monotonic increasing. 
Proof of Proposition [II' 

We exploit the permutation symmetry of the problem in the manner pioneered 
in |Pos56j . |Pos62] . and advanced in |Hal83j . |Hal88j . |Hogll| . For > 2 



the Hamiltonian ([T]) is a permutation-symmetric sum of one- and two-body 
operators, and to find the infimum of its spectrum it suffices to let it act 
on the closed subdomain of bosonic (i.e. permutation-symmetric) A^-body 
wave functions This allows one to rewrite the energy functional ([2]) 

as a conditional ^'^'^^-average of a two-body operator for all N > 2. This 
transforms the bosonic ground state problem of the Hamiltonian ([T]) into the 
realm of Helium-type ground state problems. We implement a variation of 
this theme. Using rescaled marginal and conditional probability densities of 
gi'^) = we rewrite the rescaled energy functional f|T3|) as a symmetric 

— 2-point average of a two-body energy functional which is evaluated with 
a ^AT-symmetric two-point probability density conditioned on the remaining 
N — 2 points. 

More explicitly, invoking the notation stipulated in section 12. 3[ we write 
(gi,g2) =: q^^^ e and (g3,-,giv) =: g^^~^^ G M3(7V-2)^ 
note the A^ — 2-th marginal probability density of g^'^\ i.e. Q^N-2iq^'^~'^^) '■~ 
J g^^'^{q^^^)(fq. By we denote the conditional two-point probability den- 
sity obtained from g^''\ i.e. g{^]\q^^^q^''-'^^) := Q^^Kq^"^)/ QN^iq^''''^^)- Now 
define the two-body energy functional 

This functional is of the Born-Oppenheimer Helium-type; when k = 1/2 and 
A = 1 it becomes the Born-Oppenheimer Helium energy functional rescaled 
to "normal form;" see [Thi02] . Note that by scaling q — t- X^^q one can always 

accomplish "normal form," viz. 'E(^^a,k(^'*'^'') ~ -^"^"^^x!! fc/xi^^"^^) ■ 

We are now ready to express the rescaled energy functional f|T3l) in terms 
of ([25]) and g'^''\ For A^ = 2, we obviously have £^^^^(^^^0 = '^S,x,i/2{0^^^) ^ 
while for A^ > 2 we hav^ 

= / h^'l ^ {^^\ ■ W'-'^))] 9'.^Uq'-'')d''''-l (26) 



(JV-2) 



^The redundant square brackets inside the integral are meant to facilitate the parsing of 
the formula. 
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With ( 126|) the ground state problem of the A^-body Hamiltonian ([T]) reduces 
to the discussion of averages of Born-Oppenheimer Hehum-type operators. 

Having ( |26l) . we now use that iVi— )■ 1 — 1/A^is increasing, and that the set of 
two-point conditional probability densities of S'^-symmetric A^-point densities 
is decreasing in N in the sense of set-theoretic inclusion. □ 

Remark 10. PropositionUl and its proof are essentially contained in a recent 
paper by Hogreve [Hog 1 VI , who proved for neutral Born-Oppenheimer atoms 
that E'^{N,N)/N^ grows monotonically in N > 1. Inspection of Hogreve 's 
proof, see his formulas (7)-(10), reveals that it generalizes verbatim to Born- 
Oppenheimer ions, without restrictions on X, and that our proof above (when 
X = 1) is just Hogreve's argument, recast in the probability measure setting. 

Note that the restriction > A^* was not needed in our proof of Proposi- 
tion [U it enters next. 

We now turn to the last part of Theorem [3] and show that the limit ([23]) 
exists, and that it is negative, strictly decreasing, and concave for A > A*. But 
armed with Prop. [1] this is easy. Namely, as is well-known, the ground state 
energy 8.^{XN, N) is a negative, strictly decreasing, and concave function of 
A > A* when A^ > A^* (cf. |Thi02j ). Furthermore, by our discussion of the 
upper Hartree bound on N~^8.^{XN, N) we have, uniformly in A^, 

N-'E-Z{XN, N) < inf {t^M\P ^ ^(K^) n Wp ^ ^\^')}}, (27) 

and it is easy to see that the infimum of the right-hand side is negative and 
strictly decreasing I — oo as A t oo. Pairing this with Proposition [T] gives us 

Corollary 2. The limit eoo{X) = \imf^-^ao N~^E^{XN, N) exists. Furthermore, 
£:oo('^) is negative, strictly decreasing, and concave. 

Of course, by Corollary [2] and (!27|) . also £00 (A) inherits the upper bound 
£oo(A) < inf {eoo,A(p)|p e ^(M=^) n Wp G i3^(M=^)}}. (28) 

We next show that the sequence of marginals {Qx^I} N>raa^{n n*} weak 
limit points in 9Jt"(R3") n {y/pn E i^^R^")} whenever A > A*. To this end we 
need another auxiliary result which likewise is of interest in its own right. Its 
key element is a super-additivity result by Carlen jCarQlj . 

Proposition 2. For any ^('^^ E^T] {y/g'^''^ E Sj^}, the map n ^ T^'''^ 
has the following properties: 

(A ) Positivity.- For all n< N, 

>0; (29) 
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(B) Monotonic increase; If k <n < N, then 

-^(n)(^(^))>^W(^(-)). (30) 

(C) Super-additivity; For n = k + i < N , 
Proof of Proposition [B' 

Positivity is obvious. Super-additivity was proved in jCarOl] (his Thm. 3). 
Monotonic increase now follows from super-additivity and positivity. □ 
With the help of Proposition |2] we are now able to establish 

Lemma 1. For any A > A^, and n e N the sequence {y/ f^x',!} N>ma.x{n n*} ^■^ 
weakly compact in Sj^[M.^"'). 



Proof of Lemma [II' 

Let > max{n, A^*}, then = [A^/nJ + k, where \_N/n\ is the integer part 
of N/n, so k < n. We now use super-additivity (property (C) in Propj2]) and 
then positivity (property (A) in Propj2]) of Fisher information to obtain 

> m^'^'K^i:!), (32) 

cf . the Gibbs entropy estimates in [MeSp82] , |Kie93j . (Alternately, ( 15^ follows 
by first using increase (property (B) in PropJ^J to estimate jF*-^-' (^?^^'') in terms 
of the Fisher information of the n [A^/nJ -marginal of then super-additivity 
(C) to get to the n-th marginal.) Dividing (!32|) by A^ gives 

i m ^^"H^^S) < ip^^^Hri"^) VnGN. (33) 

We now invoke the virial theorem, which for any A > A^,(A^), A^ > A^*, 
gives us 

zV^n,r)) = -£:r(-)(,r)), (34) 

and therefore 

^'M^') = -'w^':!M^')- (35) 

By Proposition [1] (also recalling our rescaling which leads to (|T3|) ). the map 
A^ I—)- N~^'Ej;^\(^g^^^) is monotonic increasing, and so, by f l5B]) . we conclude 
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that N N ^ ^F^'^\g^;^'') is monotonic decreasing. Thusjfl 

jV-i<r(-)(,(-))<^(,W), (36) 
and so, by ( 15^ . 

j^[^\!F'''K(>[:l)<M')^ (37) 
for > max{n, A^*}. But since ^ |_^J — )■ ^ as — )■ oo, we in fact have 

:FW(,g)<C(n)j(,«), (38) 

for > max{ra,A^*}. With this means that W Q^^h N>mi.^{n,N*} is 

weakly compact in the homogeneous Sobolev space i3^(]R^"). On the other 
hand, we also have that {^'An}iv>max{n Af*} ^'*(M^"), and so the sequence 
WA'',l}N>m^.{n,N*} wcakly compact in ioi(M3'^). □ 
By the familiar Sobolev embedding. Lemma [1] implies: 



Corollary 3. For any A > A* and n G N the sequence {Q^xn}N> 
weakly compact in ii^(M^"') for all 1 < p < 3n/{3n — 2). 

Clearly, by Corollary |3] the sequence A^ i— )■ gi^^ for A^ > max{n,A^*} is 
weakly compact in the set of permutation-symmetric measures 971* (M^") with 
mass < 1, and which are absolutely continuous w.r.t. Lebesgue measure. We 
will now show that any weak limit point fix^n is in fact also in *P*(M'^"), i.e. 
has mass = 1. This will follow from a characterization of the weak limit points 
in terms of an A^ = oo counterpart of the finite- A^ minimum quantum energy 
principle ( |T6l) . It will vindicate all the remaining claims in Theorem |3l 

We introduce *P'*((M^)^), the permutation-symmetric probability measures 
on the set of infinite exchangeable sequences in M^. Let {/inlneN denote the 
sequence of marginals of any /i G *^*((M^)^). The de Finetti |deF37] - Dynkin 
|Dyn53| - Hewitt-Savage |HeSa55] decomposition theorem for *P*((M^)^) states 
that for each G *P*((M^)^) there exists a unique probability measure <^(dp|/i) 
on Cp(M^), such that for each G N, 

/in= /" P^"^(dp|/i), (39) 



max{n,A''*} 



IS 



where p®"(d'g) = p(dg) ® ■•■ ® p(dg), n-fold. Note that ([21]) expresses the 
extreme point decomposition of the convex set *P''((M^)^), see |HeSa55] . If p 



^Clearly, g^^"^ is the minimizer of a Born-Oppenheimer hydrogenic ion. Thus, the kinetic 
energy per particle of the ionic ground state not only decreases with increasing numbers of 
electrons, it is in fact rigorously bounded by the explicitly computable kinetic energy of an 
"electron" in the hydrogenic ground state. 
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is absolutely continuous w.r.t. Lebesgue measure, one can identify the mea- 
sures /i, fin, and p with their densities; this notational convenience should not 
cause any confusion. We also introduce the subsets ^^^((M^)^) C *P''((M^)^), 
with p > 1, as the probability measures on (M^)^ whose decomposition mea- 
sure is concentrated on *Pp(]R^) := fl £^)(M^), which are the absolutely 
continuous (w.r.t. d'^g) probability measures on M.^ whose density is also in 
£^(R'^) (of course, the wording is redundant if p = 1). Lastly, we introduce 
the abbreviation fi{g) := Jgdfi for the /^-expectation value of g. 

For any fi G *P3((]R^)^) we now define the mean Volta energy of /i as 

:= mn !,„(- AE ^ + 1 EE (40) 

l<k<n l<k<l<n 

Clearly, by the permutation symmetry of fi, this definition yields right away 

l^,(/i) = -Xp,{\q\-')+lp2{\q-q'r), (41) 
and by the linearity of /i (p) , the presentation (l39i) now yields 

^x{f^)= [ (-Ap(|gr^) + |p^^(|g-gr^))^(dp|/i). (42) 

For any p G ^*((]R^)^) for which each y/ pn G iO"'^(R'^") we also define the 
mean Fisher information of p as limit 

J(p):= hm (43) 



Here, ir'^"^(p„), n G {1, ...}, is the Fisher information of p„, as defined in f|T8|) . 
The limit f H3|) exists or is +oo. This is a consequence of the super-additivity 
and positivity of the Fisher information; see Proposition [21 which holds also 
with G q3"(M3^) replaced by p G ^^((M=^)^). 

An analog of f H2]) for the mean Fisher information holds as well, thus 

^(p) = J J(p),(dp|p). (44) 

This is a consequence of the next proposition, which seems to be a new result. 
Proposition 3. The mean Fisher functional is affine linear. 

Proof of Proposition 0' 

We use that the Fisher functional and the Gibbs entropy functional are related 
by the heat flow equation. More precisely, let p G ^*((M^)^), and suppose 
that Pk+n{t) solves dtHk+n(t) = A(''+")pfc+„(t) for t > 0, with lim^o Pfc+n(i^) = 
Pk+n, where A(fc_|_„) denotes the Laplacian in ]R3(fc+n)_ 'jj^g^ integration over 
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the k variables in A(^k+n), using that A(fc+„) = A(„) + A(fc) and noting that 
/in = / lin+k^^^q, shows that also 9j/i„(i(:) = A(„)/i„()f:) with lim^o /in = /^n- 
Since this holds for all n G N and A; G N, we have a compatible sequence 
of equations defining G *P''((M^)^), solving dtiJi{t) = Afi{t) for t > 0, 
with lim^o /i(^) = Z^- Moreover, by the linearity of the heat flow equation, 
if fJ'^^^t) and /i^^^(t) are two heat flow evolutions in ^^((M^)^), then so is 
/i(")(t) = a^(i)(t) + (1 - for all a G [0, 1]. 

Now we use that for those /i G ^^((IR'^)^) for which ^(/i) exists, we have 

^^^"^ (/in(t)) = :r(") (/i„(t)) V t > 0. (45) 

In particular, if ^(/i*-"-') exists for a = and 1, then ^(/i*-"^) exists for all 
a G [0, 1], and then holds with /i(0) = /i^") G ^^^((R^)^). 

Next, after division by n and taking n oo, f H5l) for /i(0) = /i'-"^ yields 

A5(^(")(t)) = j(/i(-)(t)) Vt>0. (46) 

Here, 5(/i) = lim„^oo n~^S^"^ (/in) is the mean Gibbs entropy of /i, cf. |MeSp82| 
and |Kie93] : it is the adaptation to the thermodynamic mean- field limit set- 
ting of the "mean entropy of states" originally introduced in jRoRu67] for 
the conventional thermodynamic limit setting of the classical grand-canonical 
ensemble; see also jRue69] . The mean entropy of /i is well-defined whenever 
inf„gN ^"^5''"^ (/in) > C, in the sense that it might be +oo (which is the case 
iff C = oo); also C = — oo can occur (e.g., if = S^^), but will not occur with 
our Born ensemble measures. Crucially important, straightforward adaptation 
of the proof of Proposition 3 in |RoRu67j (cf. |Kie93] ) yields: 

Lemma 2. The mean Gibbs entropy functional S_[fj) is affine linear. 

Thus, 5(/i(")) = aS{fi^^^) + (1 - a)5(/iW), and so 

^5(/.H(t)) = a|5(/i«(t)) + (1 - «)A5(^(o)(t)). (47) 

Now using fj46|) with a = and 1 to rewrite r.h.s. (H71) . then equating it to 
r.h.s. fH6|) . we obtain 

j(/i(")(t)) = «^(/.«(t)) + (1 - a)^(/i(°)(t)). (48) 

Taking the limit t | proves that J is affine linear for all /x^") G q?*((M^)^) 
for which ^(/i^°-') and j(/i^^-') exist. Proposition 1 is proved. □ 
Although we won't need it for our proof, the listing of properties of mean 
Fisher information (j33]) would be incomplete without the next Proposition, 
which is proved by adaptation of a proof in |RoRu67] . their proposition 4. 

Proposition 4. The functional fx h-> ^(/i) is weakly lower semi-continuous. 
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Finally we define the mean quantum energy of fi & *Pp((M^)^) by 

£oo,A(/"):=|r^(/^) + ^'eVA(/^)- (49) 
By (jS]) and (02]) we have 

^oo,a(/^) = / eoo,A(p)?(dp|/i), (50) 

J3) 



where eoo,A(p) is the quantum energy functional introduced in (l22l) . It is well- 
defined for those p e ^{R^) for which G Sj^{R^), for all A > 0. This 
definition extends by scaling to all p G 9Jl(]R'^) for which -^/p G i3^(]R^), for all 
A > 0. In all other situations it is defined as Coo^x{p) = oo. This completes our 
setup, and we are ready to characterize the weak limit points of the sequences 
of marginal Born ensemble measures. 

In the following we again invoke the restriction that > max{n, A^*} for 
each A > A* and each n G N. 

By Lemma [1] and its Corollary [3l we can extract a subsequence, denoted 
{ox^'^'n ^ ^QS^'fM^"), for which Iv/fPiTi^''} , converges weakly in 

L.CA,n J JV>max{n,]V*} ^ ^ L V t^A.n J ]V>max{n,]V* } ° 

i^^(K3"), SO that {^^r"}^>^,,{„jv.} ^^^^^^ converges weakly in ^^(M^")^ for all 
1 < p < 3n/(3n - 2). We seT 

lim = f^in- (51) 

As remarked earlier, at this point we only know that y/p\n & SS^iR^^) and 
that G aJt"(M3"), with ^^^^(ffi^") = Aa,i(^^) < 1- However, we can define 
PA,n := p\j p\^AM?); clearly, AA,n(K'") = l- 

We now use (l33l) . with A^[A^] in place for A^, and let A^ — oo. By the lower 
semi-continuity of the Fisher information |BBL81] . 



while i 



liminf:rW(,fj^l))>iFW(Al,J, (52) 
-^^ -, so from fl33|) we obtain 



N 



lim i:r(^)(,r)>ir('^)(AU' VnGN. (53) 



Here we used that we already showed that A^~^ir^^-* (^'a^'') has a limit; if that 
were not yet known, we could simply replace the "lim" by "liminf" at l.h.s.f l53p . 
Next, permutation symmetry yields the identity 

M"^(iEE k,-%\-') = {i-j,)^,st{iY.Y. (54) 

^ l<J<fc<JV ' ^ l<l<fc<7l ' 
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Weak lower semi-continuity of the Coulomb kernel, and its positivity together 
with n — 1 < n, plus 1/N ^ 0, altogether now yield, for each n G N, 

li^i^f ^f^rf^EE \Qj-%\-')>tiAMJ:T. \Qj-Qk\-')- (55) 

N^OO ^ l<j<k<N ' ^ l<j<fc<n ' 

Again by permutation symmetry, we also have the identity 

i^'(-AE.^)=i.£(-AE (56) 

^ l<fe<JV ' ^ l<k<n ' 

This time, weak 2F convergence with p = 3 yields 

^.m ir^f'(- AE = h)- (57) 

W— >00 ^ l<fc<7i ' ^ l<k<n ' 

Altogether, (jSS]), ([55]), (jSTD yield, for each n G N, that 

lim3<U^^r^)>^<\Kn) (58) 



for any weak limit point Here we used that we already showed that 

^~^'^oo!a(^'a^^) a limit, namely e^{X)\ otherwise, we could simply replace 
iV by iV[iV] and "lim" by "liminf" at l.h.s.([58D. As for r.h.s.([58l), scaling yields 

<UaU=A1i(»')<UAaJ- (59) 
Thus, for any subsequence N\N\ such that (1^ holds, we have 

e^{\) > A^,i(M3)^<\(AA,n) V n G N. (60) 

Next, {fix^n\n £ N}, with each y/fix^n ^ , is a compatible sequence of 
marginal probability measures which defines a fi\ & *P'^((R'^)^). Thus, taking 
the supremum over n of (1601) . which is equivalent to taking n — )■ oo, recalling 
the definition (H9|) of x extremal decomposition (!50|) . we obtain 



£oo(A) > Al,i(M^) y eoo,A(p) ^(dpI^A), (61) 

for any subsequence A^[A^] such that ( 15T|) holds. But 

y c^M^idplfix) >mf{t^M\p e^;^p e Sj'}, (62) 

so we finally obtain 

£oo(A) > A^,i(M=')inf{eoo,A(p)|pG^;v/pei3^}. (63) 
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We now put everything together. First of all, the estimates ( 163|1 and (128|) 
together imply that 



£oo(A) >A1,i(M3)£^(A), (64) 
but since eoo{X) < by Corollary [2l (IMI) implies that fi\ i(IR^) = 1, and so 

lim gi"^ = ilx,n V n G M. (65) 
Having /i^ i(IR^) = 1, the estimates (jHSD and (125]) together furthermore imply 

£oo(A) = inf{eoo,A(p)|p e q}; G -^i'}- (66) 
But then, by fl6Tl) and fl62l) . we also have 

y eoo,A(p)?(dp|AA) = inf{eoo,A(p)|pe^;v/pef|i} (67) 

for any weak limit point fi\ of the rescaled sequence of Born's ground state 
ensemble measures. It follows that the decomposition measure (;{dp\fix) is 
supported on those p G ^ which actually minimize eoo,A(p); for suppose not, 
then the average / eoo,A(p) ^{dplfj^x) > inf{eoo,A(p)|p G a/p G Sj^}, in contra- 
diction to (ETj). Therefore the infimum is in fact attained; furthermore, since 
A//i„ G i^^, we can conclude that the minimizer is actually in^n{^pGi?^}. 

Lastly, by the well-known convexity of p i— )■ !F{p) (see |BBL81j ). the mini- 
mizer is unique. By the weak compactness of the sequences of marginals this 
in turn implies that all subsequences converge to the same limit. 

This proves Theorem [31 □ 
This proves Theorem [H □ 



3.2 Proofs of Theorem [2] 

Theorem [2] becomes a corollary to the proof of Theorem [HI once we add the 
random variables into the story. This step supplies proper meaning to the 
notion of l?/'*^^'' p as density of a joint probability measure of an N-body system, 
which otherwise would be just another name for a positive function which 
integrates to unity. As Born emphasized, it is this step which relates ■ip^'^^ 
to the empirical "dots on the screen" gathered in laboratory experiments. An 
even more profound insight into the physics is obtained by thinking of Iz/''-^'* ^ as 
density of a typicality measure. In this section we pursue both interpretations 
of 'ip^^\ Also here we have to impose that > A^*. 



22 



3.2.1 Probability 



So we finally enter the probabilistic meaning of Born's ensemble for the rescaled 
ground state if)'^^^ of the Hamiltonian ([T]), namely family {Q^^\k e Nj 
of i.i.d. copies of a random vector Q^^-* e with normalized, stationarjo 
ensemble probability density q"^^ = l'?/'^^^^. Physically it is more natural to 
think of the vector Q^^^ G M'^^ as a random set of points {Qi, Qjv) 

1. e. as a representative of the equivalence class of such vectors under the 
permutation group Sn, also denoted by Q^^\ Associated to each random set 
of points Q^^^ is a unique normalized empirical one-point random "density" 

4'(-)(^) = ^ E SqM), (68) 

l<j<N 

and more generally a normalized empirical random [/-statistic ^q^n) defined 
in ( JT2I) . The quantities of interest to physics are of the form 

Prob(rfKR(A(;i),p(")) ><5), (69) 

where "Prob" refers to l^/'^^^P as probability density for Q*-^-* G M'^^, and 
p(")(si, Sn) G (^P** n (r°)(M^"') is some continuum approximation to a phys- 
ically meaningful macrostate (see below); in practical situations, n = 1 or 

2. Furthermore, c/kr is some Kantorovich- Rubinstein metric quantifying the 
weak topology on *p*(M^"). 

We inject that we have the luxury of being allowed to remain somewhat 
vague about both, the notion of macrostate in general, and which Kantorovich- 
Rubinstein metric we want to use. Namely, since we are invoking the physi- 
cists' notion of "ideal position measurements," it is immaterial how approxi- 
mately to these one chooses a "macrostate" in real life, and we can allow any 
(piecewise) continuous ?7,-point function as "macrostate;" for more on conven- 
tional macrostates, see |GoLe04] . And since the A^-body ground state wave 
function decays exponentially in R^'^ Agm82| , and since the same is true in 



for the Hartree minimizer |BBL81] , when paired with our Theorem [T] we 
can conclude that the exponential decay is uniform in A^ for any nth marginal, 
so that the particular choice amongst the standard cost functions in the KR 
metric is of no essential importance. 

For any fixed n < N, with A^ > A^*, Born's ensemble entails the conven- 
tional law of large numbers. Namely, let !N be the number of ideal measure- 
ments of Q^^^ that are being performed, with A^ > A^* fixed. Then the law of 
large numbers says that when Jsf — )■ 00 the sample mean over the independent 
and indentically distributed (i.i.d.) empirical A^-point densities converges in 



^The density is stationary in the sense that the wave function -0^"-' g£^(R'^^) is 

stationary under the action of the unitary group generated by the rescaled Hamihonian ([l} . 
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probability to the theoretical mean, given by IV'a so, for any 

n < N, with N > N*, the sample mean over the i.i.d. n-point densities 
converges in probability to their theoretical mean 

ri^fe,-,0= / l^f¥(g("^)dg„+i---dg^. (70) 

Our Theorem [T] then adds to this that the theoretical mean d'xni^i^ ■■■^^n) 
takes a simple Hartree product form when N ^ oo. 

Our Theorem [2] says something much stronger, namely that, in probability, 
the outcome of a single ideal measurement of the empirical n-point density of 
an A^-body system will converge to the Hartree product density when N ^ oo. 
The number !N of ideal measurements of the empirical n-point density plays 
no role in the argument. Of course, since a single ideal measurement of the 
empirical n-point density of an A^-body system is equivalent to ideal mea- 
surements of the individual particle positions, also Theorem [2] can be phrased 
as a law of large numbers, this time the number of particles in a single 
system. 

A first proof of Theorem [H" 

We recall that, if Q^^^ denotes an infinite random sequence of points in 
which are independently and identically distributed by a single-point prob- 
ability measure p (having density p), then their joint probability density is 
the infinite product p*^^. Considering Q^^^ as a subsequence of Q^^\ the con- 
ventional law of large numbers states that the associated sequence of n-point 
empirical measures ^q(n) converges in probability for each n G N; viz. for 

each n G N: c^kr ( Aqcat) ? P®") 0, in probability. 

Now, as explained in the proof of Theorem [3], the N ^ oo limit of Born's 
ensemble measure yields g^^^ pf^, an infinite product measure on the set 
of infinite sequences q^^\ 

This proves Theorem |21 □ 
The proof of Theorem [2] just given brings the law-of-large-numbers charac- 
ter of Theorem [2] to the fore, but leaves the interpretation in terms of "a single 
measurement of the empirical n-point density of an A^ = oo system" rather 
implicit. To bring out this aspect of Theorem |2l we now give a second proof. 

A second proof of Theorem\^ 

We already emphasized that each q'^^^ C M'^ is associated with a family of 
normalized empirical n-point densities {A^"lr)}n<]v Now consider g^^^ as a 
finite subsequence of an infinite sequence g^^^ of points in M'^. Each infinite se- 
quence g*^^^ C is associated with a family of sequences {A^"lr)}jv>max{n,iv*} £ 
^''(R^"), n G N. For each g^^^ for which there is a p G "^{M?) such that 
A^{iv) P weakly in the sense of measures, we can define ^^il) to mean the 
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normalized measure p G ^(M^) given by the weak limit of ^Jn) when N oo. 

Each such A^(],) = p is just the first member of an infinite family of normalized 

n-point measures A^"i), G N, which are defined analogously as weak limits 

of the n-point empirical measures A^"]r), G N, as — t- oo. It is easy to see 

from the definition of A^"]r) that, inevitably, A^"i) = p®", n G N. 

Now, for each n < N, with N > N*, Born's ensemble measure \ip{^''\'^d^q 
can be identified with a probability measure <;-|"'^'' on the set of normalized n- 
point measures; note that lives on the singular subset given by the generic 
version of 1^, viz. ([I2]) with Q^'^) replaced by q^""^ = G M^", or 

rather {q^, g^} C M^. By the explanations given in the previous paragraph, 
whenever (^^^-^^ has a limit as — )■ oo, then <j^"'' is a probability measure on 
the set of normalized n-point measures which lives on the subset of these which 
are n-fold product measures, viz. <^|"^ is supported by the set of p®"". Clearly, 
all the <;^\ n G N, are uniquely determined by <;}^^ =: ^a- Note that this is 
precisely the content the Hewitt-Savage decomposition theorem! More to the 
point, the proof of Theorem [3] shows that the Hewitt-Savage decomposition 
measure ^(dplp^) simply is the limit — )■ oo of Born's ensemble measure 

|^^^v)|2^3Ar j^gj^^jgg^ with ^l^'^"*. 

Thus, the fact, found in the proof of Theorem [31 that <j(dp|pA) is concen- 
trated on a single px (which is the unique minimizer of the asymptotic Hartree 
functional), means that the probability, of picking, from Born's ensemble, an 
empirical configuration which differs from p\, goes to zero when — > oo. 

This is the content of Theorem [2l □ 



3.2.2 Typicality 

Born's statistical / probabilistic interpretation of {ipxl"^ inevitably suggests the 
interpretation of Theorem [2] in terms of a law of large numbers. This notion 
refers to the sample mean over many experimental results becoming sharply 
distributed. Yet, in our second proof of Theorem [2] we already saw that we 
can make sharp statements about the outcome of a single measurement of the 
n-point densities in individual systems with overwhelming likelihood. 

Even more penetrating is the notion of typicality, in which IV^^^"*!^ features 
as having the meaning of tj^icality measure rather than probability measure. 
The notion of typicality also refers to individual systems, but it refers to these 
in their own right, without reference to any measurements being performed 
on them: If the overwhelming amount of configurations Q^^^ correspond to 
the same p|yp, then this p^yp is characteristic, i.e. typical for the A^-body 
system. It is this notion which really lies at the heart of Boltzmann's insights 
|Bol96] about the ergodic ensemble, not the subsequently invented dynamical 
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notion of "ergodicity;" see |Gol01] . and also |Kiell] . A typical n point density 
Ptyp(si, Sn) G (^"^ n is therefore defined implicitly as the — in the 

simplest case: unique — function for which 

Meas {dKK ( aJ),, , pi^^) > 5) V<5 > 0, (71) 

where "Meas" refers to |V'a^''P as typicality measure (density) on M.^^ . 

Clearly, everything we said about Theorem |2] in terms of probability can 
be rephrased in terms of typicality. □ 



4 Concluding remarks 

In this paper we have presented a novel approach to the Hartree limit of 
bosonic Born-Oppenheimer atoms or ions. Some of the items listed in our 
Theorem [1] have previously been proved with different techniques, while other 
items seem new, indeed. Our approach is based on Born's statistical ensemble 
interpretation of see Theorems [2] and [31 Their proofs are inspired by an 
established strategy in classical statistical mechanics, pioneered in |MeSp82] . 
The key idea is that the role played by the negative Gibbs entropy in the clas- 
sical petit-canonical ensemble, is played by the Fisher information functional 
in Born's quantum ensemble. 

Along the way we also established an auxiliary result, namely the mono- 
tonic increase of N~^E'^{XN, N) with A^, see Prop. [H The monotonicity 
of N~^S,'^{N, N) (i.e. A = 1) had already been shown by Hogreve Hogll , 
who used it to establish the existence of a stable bosonic "He~ ion" in Born- 
Oppenheimer approximation. 

Due to our currently incomplete knowledge of the optimal range of Z val- 
ues for which a proper bosonic ground state "^[^'^^ exists, the proofs of our 
theorems work under the restriction that A^ > A^*, where A^* is some finite 
positive integer. In this regard we formulated two reasonable conjectures, 1 
and 3, which, if true, would imply that A^* = 1, thereby eliminating the need 
to restrict N to N > N*. Our conjectures are also of some interest in their 
own right; in particular, as pointed out by the referee. Conjecture 1 is closely 
related to a known conjecture, 2, formulated in [LiSelOj . It would be very nice 
if all these conjectures would be settled in some future work. 

As also noted by the referee, the convergence of the states accomplished 
here is relatively weak, and it would be very interesting to know if the technique 
could be modified to prove the convergence in trace norm of the family of 
reduced density matrices for the ground states. 

Other projects worthy of our attention in the future include the general- 
ization of our technique to the conditional ground state ensemble of Galilei- 
invariant "bosonic atoms and ions;" i.e. to drop the Born-Oppenheimer ap- 
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proximation and to treat the system in its center-of-mass frame. In the same 
category are "bosonic neutron stars." Less obvious, and therefore more in- 
triguing, is the generahzation to fermionic atoms, ions, and stars. In this case 
one would expect to approach Hartree-Fock theory |LiSi77] . 

In this work the terminology "bosonic atom" invariably meant an atom 
whose "electrons" were bosons, and an atom with true electrons was called 
"fermionic atom." Same for ions. While this terminology has been in use by 
parts of the physics community for some time by now, the term "bosonic atom" 
can also be found in the physics literature to mean an atom with conventional 
electrons, but whose total spin angular momentum quantum number is an 
integer. The most prominent example, perhaps, is the neutral helium atom 
^He in its ground state, having spin 0. Experimental studies of ^ 1 such 
bosonic atoms in a trap have been making headlines in the recent past with 
the empirical observation of Bose-Einstein condensation in gases consisting of 
bosonic alkali atoms |Aetal95t IDetal95] . 

Since theoretical physicists have employed the Hartree approximation also 
to study Bose-Einstein condensation, it is worthwhile to point out that the 
mathematical approach to the Hartree limit developed in the previous sections 
readily extends to the Hamiltonian used by theorists studying BE condensates, 
which yields the energy functional ( 1T3|) with ze = 1 and with potential energy 
functional ( IT5l) replaced by 

= /(a E + ^ E E U{\q, - q,\)) ,^-^d\ (72) 

J ^ l<fe<JV l<k<l<N ' 

here, A > now represents the strength of the confining harmonic potential, 
and [/ > is some short range, integrable repulsive potential (frequently even 
taken to be bounded). We leave it to the interested reader to verify that our 
Hartree limit theorems and proofs extend nearly verbatim to this setting; some 
of the technical estimates even simplify. 
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